Accepted for publication in ApJ. 

Preprint typeset using I^T^^X style emulatcapj v. 8/13/10 



CONSTRAINING DUST AND MOLECULAR GAS PROPERTIES IN Lya BLOBS AT z - 3 

YuJiN Yang\ Roberto Decarli^, Helmut Dannerbauer^'^, Fabian Walter^ Axel Weiss'", Christian Leipski\ 
Arjun Dey^, Scott C. Chapman'', Emeric Le Floc'h^, Moire K. M. Prescott'', Roberto Neri**, Colin Borys'^ 
YuiCHi Matsuda"\ Toru Yamada"\ Tomoki Hayashino^^, Christian Tapken^^, Karl M. Menten"* 

Accepted for publication in ApJ. 

ABSTRACT 

In order to constrain the bolometric luminosities, dust properties and molecular gas content of giant 
Lya nebulae, the so-called Lya blobs, we have carried out a study of dust continuum and CO line 
emission in two well-studied representatives of this population at z 3: a Lya blob discovered by 
its strong Spitzer MIPS 24^m detection (LABd05; Dey et al. 2005) and the Steidel blob 1 (SSA22- 
LABOl; Steidel ct al. 2000). Wc find that the spectral energy distribution of LABd05 is well described 
by an AGN-starburst composite template with LpiR = (4.0±0.5) x 10^^ L©, comparable to high- 
z sub-millimeter galaxies and ultraluminous infrared galaxies. New Large APEX Bolometer Camera 
(LABOCA) 870/.tm measurements rule out the reported Submillimcter Common-User Bolometer Array 
(SCUBA) detection of the SSA22-LAB01 (5*850^01 = 16.8 mJy) at the > 4cr level. Consistent with this, 
ultra-deep Plateau dc Bure Interferometer (PdBI) observations with ^2" spatial resolution also fail to 
detect any 1.2mm continuum source down to « 0.45 mJy bcam~^ (S"")- Combined with the existing 
(sub)mm observations in the literature, we conclude that the FIR luminosity of SSA22-LAB01 remains 
uncertain. No CO line is detected in either case down to integrated flux limits of S^/SV < 0.25 - 1.0 
Jy kms~^, indicating a modest molecular gas reservoir, M(H2) < 1-3 x IQ^^Mq. The non-detections 
exclude, with high significance (12 cr), the previous tentative detection of a CO J=4-3 line in the 
SSA22-LAB01. The increased sensitivity afforded by the Atacama Large Millimcter/submillimcter 
Array will be critical in studying molecular gas and dust in these interesting systems. 
Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium — radio 
lines: galaxies — submillimcter: galaxies 



1. INTRODUCTION 

Lya nebulae (or Lya "blobs" ) are extended sources at 
z 2 - 6 with typical sizes of the Lya emission region 
of ^ 5" (> 50 kpc) and line luminosities of LLya ^ 10'''^ 
ergs~^, and are so me of the most mysterious astr onom- 
ical objects ( e.g.. iKeel et al.l Il999t iSteidel et al.l 
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emission is poorly understood. For example, Lya blobs 
may represent galaxies forming via cold gas accretion 
dFardal ct al. 2 001: Haiman. Siiaans. & Ouataertl 120001 : 
iDiiicstra fc"Loebll2009l:lGoerdt et al.ll2010i). galactic-scale 



outflo ws due to star formation (jTaniguchi fc Shioval 
|2000[) . or the result of intense r adiative feedback from 
active galactic n uclei (AGN) (jHaiman k. Reed 120011 : 
IGeach et al.1[2009t) . 

To understand the nature of these Lya blobs, we first 
need to constrain their energy budget, i.e., the bolomet- 
ric luminosities of the galaxies within or in the vicin- 
ity of the Lya halos and compare the available energy 
with the observed Lya luminosities. Many Lya blobs 
do not appear to contain bright UV continuum sources, 
bringing into question the nature of the underlying power 
source for these nebulae. Infrared (IR) and (sub) mm ob- 
servations provide complementary information in such 
cases, by detecting and constraining the luminosity of 
any dust-enshrouded and obscured power-sources (star- 
forming or AGN) within the nebula. So far, a handful 
of Lya blobs have been detected in the far-IR (FIR), 
suggesting that at least some Lya blobs contain energy 
sources that can power the entire Lya luminosity if only 
a few percent of their bolomet ric lumin osities ar c con- 
verted to Lya radiation (Chapman "eFaI|[20 04: Dc v^eFaLl 
[20051: IGeach et al.l[2005lr2009HColbert et al.ii200a ifnTT 
On the other hand, it appears that some Lya blobs do 
not contain any obvious energy sources detectable at 
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mid-IR or millim eter wavelengths (jNilsson et alj 120061 : 
iSmith et al.l 120081 ): in these cases the Lya emission has 
been attri buted to gravitational cooling due to cold m ode 
accretion (|Keres et alJl2005L [2009t lOekel et al.llMl . 

Despite the enormous star formation rates indicated 
by the FIR luminosities of some Lya blobs (up to ~ 
lOOOM0yr~^), their molecular gas content - where the 
stars should form - has been mostly unconstrained. 
Therefore, it is important to understand whether there 
is enough cold molecular gas to power star forma- 
tion in Lya blobs or whether the molecular gas is de- 
pleted rapidly due to galactic-scale feedback either by 
superwind-driven shocks or photoionization by AGN. In 
addition, if detected, carbon monoxide (CO) lines that 
trace the molecular gas can be used to probe the kine- 
matics and the excitation conditions of the surrounding 
gas. Lya blobs also present an unique opportunity to in- 
vestigate the molecular gas content in environments that 
differ significantly from sub-millimeter galaxies (SMG), 
high redshift quasars (QSOs) or ultraluminous infrared 
galax ies (ULIRGs) that have been s tudied in the past 
(e.g., iSolomon fc Vanden BoutI I2005D . Most luminous 
Lya blobs share the common pr operty that they re- 
side in overdense e n vironments (iPalunas et al.l 120041: 
Matsuda etlUI 120041: [Prescott et a l.l 120081 lY ang et al.l 
20091 120101 iMatsuda et al.l 120091 ). suggesting that they 
may represent the sites of massive galaxy formation. Fur- 
thermore, the discovery of spatially-extended He II A1640 
emission and very weak metal emission lines from a Lya 
blob at z = 1.67 sug gests that the gas in the blob may 
be of low metallicity (jPrescott et al.ll2009( ) . 

In this paper, we present new (sub) mm obser v ations 
of a z = 2.656 Lya blob discovered bv lDev et al.l ()2005l ) 
and measure its bolometric luminosity and dust prop- 
erties. This blob (SST24 J1434110 +331733; hereafter 
LABd05) was discovered i n the NOAO Deep an d Wide 
Field Survey Bootes field (jJannuzi fc Devlll999t) by its 
strong Spitzer Multiband Infrared Photometer (MIPS) 
24/.im flux, indicating that as an IR-bright source it may 
be an ULIRG with a bolometric luminosity as large as 1- 
8 X IO^^Lq. While a 350/^m flux density of LABdOS had 
been measured as a part of a large Caltech Submillime- 
ter Observatory (CSO) Submillimeter H igh Angular Res- 
olutio n Camera (SHARC-II) program (jBussmann et al.l 
[200l . the lack of (sub) mm observations covering both 
sides of the thermal dust peak of its spectral energy dis- 
tribution (SED) prevented a direct constraint on its bolo- 
metric luminosity. 

The best-studied Lya blob (SSA22-LAB01 at z = 
3.09) in the SSA2 2 protocluster region discovered by 
iSteidel et al.l ((20001) has been reported to be associated 
with a bright SMG from the Submillimeter Common- 
User Bolometer Arr ay (SCUBA) observations (e.g., 
iChapman et al.ll2004[ ). However, follow-up observation 
with the Submillimeter Array (SMA) at higher spatial 
resolution (^^2") yielded a non-detection (jMatsuda et al.l 
|2007| ). Furthermore, an Atacama Submillimeter Tele- 
scope Experiment (ASTE) AzTEC single-dish observa- 
tion with a larger beam size (w 28") than SCUBA (15") 
also failed to detect this source down to ^ 3 mJy (3cr) at 
1.1mm (|Kohno et al.ll200l iTmlura et al.l [20091 ). Thus, 
the FIR luminosity of SSA22-LAB01 and the exact lo- 
cation of the dust continuum within the blob is un- 



certain. In this paper, we revisit SSA22-LAB01 with 
Large APEX Bolometer Camera (LABOCA) measure- 
ments whose wavelength and beam size are close to those 
of the previous SCUBA observations. We also present 
ultra-deep mm intcrfcromctric continuum observations 
aiming to detect the possible compact sources within the 
blobs and, thus, to pinpoint the exact location of the 
energy source. 

Prior to this work, the only blob targeted wi t h CO 
observations was SSA22-LAB01. IChapman et"an (pOOl 
reported a tentative detection of a CO(4-3) line, which 
implied a significant molecular gas reservoir, M(H2) ~ 
lO^^AfoQln this paper, we present a sensitive search for 
CO emission in this and the Lya blob, LABd05. In fj2l 
we describe our dust and CO observations. In H3.ll we 
present the SED of LABd05 and constrain its dust prop- 
erties. In H3.2\ we present the singlc-dish submm obser- 
vation to verify the previous FIR measurements, and the 
deep mm observation of SSA22-LAB01 conducted to de- 
termine the location of its submm continuum emission. 
In §3.3l and §3.41 we put constraints on the molecular gas 
contents of both Lya blobs and compare their FIR and 
CO luminosities with those of other high-z galaxies. In 
Sj4j we summarize the results. Throughout this paper, 
we adopt the following cosmological parameters: Hq = 
70kms-i Mpc-\ ilu = 0.3, and = 0.7. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Observations of LABdOS 

2.1.1. IRAM-SOm 1.2mm Observation 

We observed LABd05 with Max-Planck Milhmeter 
Bolometer (MAMBO-2), on the Institut de Radio As- 
tronomic Millimctriquc (IRAM) 30 m telescope in 2005 
May 18, 19 and 22. MAMBO-2 has a half-power spectral 
bandwidth from 210 to 290 GHz, with an effective band 
center of ^250 GHz (1.20 mm) and a beam full-width at 
half maximum (FWHM) of 10'.'7. LABd05 was observed 
in on-off observing mode under good weather conditions 
for ~3 hr down to an rms of only 0.34 mJy (1 a) and was 
detected with a flux of 2.66 mJy (i.e., at 7.8 cr). 

2.1.2. New Spitzer MIPS 70^m and 160^m Observations 

iDev et all ()2005l ) reported non-detections of LABdOS 
at 70/xm and 160^m based on the shallow observations 
of the Bootes field carried out with the MIPS instrument 
on-board Spitzer. To further constrain the photome- 
try of the source at these wavelengths, deeper pointed 
observations at 70/im and 160/xm were obtained with 
the MIPS "photometry mode" as part of a General Ob- 
server program (ID 20303). The point spread function 
sizes (FWHM) are 18" and 40" at lOfim and 160/^m, 
respectively. We reduced these new observations with 
the MIPS Data An alysis Tool foll owing the prescriptions 
described by Gord on et al.l ()2007t ) , and derived fiuxes for 
L ABd05 using aperture photometry at the position of the 
object. However, no detection was obtained even in the 
deeper data. We estimated the 3a upper limits as the 
dispersion of fiuxes measured in similar apertures ran- 

^ Wc converted the reported CO flux to the molecular gas mass 
using the same as sump tions, cosmological parameters adopted in 
this paper fSection l3.3| l a nd the CO(4— 3)/CO(l— ) brightness ratio 
of ~ 0.5, i.e., the same as IChapman eFaLl l(200l ). 
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Hcrschel/SPIRE 250/.tm, 350/.tm, and images of LABdOS. The location of the MIPS 24^m source is marked with the 



domly placed over the background of each image. Aper- 
ture corrections were applied following the MIPS Data 
Handbook of the Spitzer Science Center, which led to 3cr 
limits of 9 mJy and 51 mJy at 70/um and 160/xm, respec- 
tively. 

2.1.3. Archival Herschel SPIRE BSO/im and SSOfim data 

The Herschel Spectral and Photometric Imaging Re- 
ceiver (SPIRE) data were obtained in SPIRE/PACS par- 
allel mode as part of the Herschel Multi-tiered Extra- 
galactic Survey (HerMES) key project (PI: S. Oliver). 
Among the observations available for the Bootes field, 
a total of seven individual maps cover the source posi- 
tion. The data were processed in a standard fashion us- 
ing Herschel Interactive Processing Environment (HIPE) 
7.0 and the latest calibration files. While the seven obser- 
vations were processed and calibrated individually, the 
resulting final map was created by combining all cali- 
brated data and mapping all the scans simultaneously. 
We then used a S PIRE source extractor implemented 
in HIPE (based on lSavage fc O lived [2007[) to find unre- 
solved sources and measure their fluxes, using data for 
the approximate shape of the beam from the SPIRE Ob- 
servers Manual (version 2.3). Photometric uncertainties 
were determined by measuring the pixel-to-pixel rms in 
an area of 180 arcmin^ around the science target. This 
rms was determined in a cleaned map, which had all the 
sources detected in the previous step removed. We show 
the SPIRE 250, 350, 500/im images in Figure □ The 
FWHMs of beam profiles are 18", 25", and 36", respec- 
tively. In the 350/xm and 500/im bands, the LABd05 
is blended with the neighboring galaxy at the south- 
west direction. It was not possible to reliably obtain 
the 500/im photometry due to severe blending. Our 
350//m SPIRE photometry (27 ± 5 mJy) is consistent 
with the previo us CSO/SHARC-II ob servation (S'ssopm 
= 37± 13 mJv: iBussmann et al.|[2009| ) within the uncer- 
tainties. We summarize the photometric observations in 
the Table m 

2.1.4. IRAM-SOm CO Observations 

LABd05 was observed with the single pixel heterodyne 
receiver on the IRAM-30m telescope in UT 2005 May 18 
and 21 in good weather conditions. We used the AB and 
CD receiver setups, with the AB receivers tuned to CO 
J=3-2 (3 mm band, Vohs = 94.614 GHz) and the CD re- 



ceivers tuned to CO J=5-4 (2 mm band, I'obs = 157.675 
GHz). The FWHM of beam is 26" and 16", respectively. 
Average temperatures were 195 K at 2 mm and 135 K at 
3 mm. Data were taken with a wobbler rate of 0.5 Hz 
and a wobbler throw of 50" in azimuth. The pointing 
was checked frequently and was found to be stable within 
3" during the runs. Calibration was done every 12 min 
with standard hot/cold load absorbers, and we estimate 
fluxes to be accurate to within ^-^10-15% in both bands. 
We used the 512 x 1 MHz filter banks for the 3 mm re- 
ceiver and the 256 x 4 MHz filter banks for the 2 mm 
receivers. As part of the data reduction, we dropped all 
scans with distorted baselines, subtracted linear base- 
lines from the remaining spectra, and then rebinned to 
a velocity resolution of 40kms^^. The total on-source 
time is 3.4 hr for both CO observations. The conver- 
sion factors from K (Ta* scale) to Jy at our observed 
frequencies are 7.0 JyK~^ and 6.1 JyK~^ for 2 mm and 
3 mm bands, respectively. The resulting rms noises are 
0.4mK (2.8 mJy) and 0.45mK (2.8 mJy) for 2 mm and 
3 mm bands observations, respectively. We summarize 
the CO observations in the Table [51 

2.2. Observations of SSA22-LAB01 
2.2.1. APEX LABOCA 870fim Continuum Observations 

The dust continuum of SSA22-LAB01 was ob- 
served at 870Mn i using Large APEX Bolometer Cam- 
era (LABOCA; iSiringo etall 120091 ) at the Atacama 
Pathfinder Experiment (APEX) in 2011 August 18. Ob- 
servations were carried out in the photometric on-off 
mode under good observing conditions (PWV — 0.5- 
0.7 mm). Throw and speed of the chopping secondary 
were set to 60" and 1.5 Hz, respectively. Chopping was 
carried out in symmetric mode with a nodding time of 
20 sec. LABOCA's spatial resolution is 19"2, comparable 
to the resolution of pr evious observations with SCUBA 
(j Chapman et al.l l2004f ). Calibration was achieved us- 
ing Uranus (66.2 Jy at the time of observations) and 
the absolute calibration was estimated to be about 15%. 
Pointing was checked every 40min on the nearby QSO 
J2225-045 and found to be stable within 3'.'6 (rms). 
The data were reduced using the BoA software package. 
The total integration time on SSA22-LAB01 is 2h20min 
(on+off, including instrumental overheads). Our obser- 
vations yield a non-detection of the source with Ssrofim 
= -2.9 ±2.7 (±4.0) mJy beam"!. The first error is the 
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Figure 2. Spectral energy distribution of LABdOS from the 
rest-frame UV to radi o . Th e squares represent the photometry 
compiled by IDev et ahl ||2005I 1 and the large circles arc the Spitzcr 
MIPS 70fj.m, leOAtm, Herschel SPIRE 250^lm, 350/i.m, and IRAM- 
30m/MAMBO-2 1.2 mm photometry from this work. The error 
bars and upper limits are given at 3a. The solid and dot-dashed 
hues represe nt the SEP templa tes of IRAS12112 -1-0305 (starburst- 
dominated: IRieke et al.l 120091) and Mrk 231 (AGN-dominated; 
IPoUetta et al.l 120071) . Templates are scaled to match the 3.6/xm, 
4.5/j m, and S.S/jm f l ux den sities. The inset shows the IRS spectrum 
from lColbert et al.l 11201 II V The AGN template shows good agree- 
ment with data while the starburst template (IRAS12112 -1-0305) 
over-estimates the FIR flux by a factor of '^5. 

rms calculated from error propagation via the rms on the 
reduced bolometer time line. The second error in paren- 
thesis is the uncertainty calculated from the dispersion 
of the individual scans, therefore represents a more con- 
servative error estimate. 

2.2.2. PdBI 1.2mm Continuum Observations 

We observed SSA22-LAB01 between May 2010 and 
January 2011 with the Plateau de Burc Interferometer 
(PdBI) in D configuration (Project ID: TGB2). The 
on-source observing time, corresponding to the full ar- 
ray of six antennae, was 5.9 hours. We set the phase 
center at the reported SCUBA submm galaxy (SMM 
J221725.97+001238.9). The 1mm receiver WideX was 
tuned at 241.5 GHz, corresponding to 1.25 mm in the 
observed and ~400 /im in the rcst-framc. The total 
bandwidth of our dual polarization mode observations 
was 3.6 GHz. The data were calibrated through ob- 
servations of standard bandpass (3C 454.3, 2223-052), 
phase/amplitude (2131-021, 2145-h067, 2223-052) and 
flux calibrators (MWC 349), and reduced with the 
GILDAS software packagesQ CLIC and MAPPING. The final 
map, created using natural weighting, has an rms noise of 
0.15 mJy integrated over the full bandwidth. Note that 
this new observation is approximately four times deeper 
than th e previous S MA observations ((Jrms = 1.4 mJy at 
880^m: lMatsuda e t al.. 20071 and one of the deepest ob- 
servations ever carried out at ~lmm on a single target. 
The FWHM of the beam is 2.4"xl.6" (18 kpcx 12 kpc) 
at 1.25mm (P.A. = 14.9°), similar to that of the SMA 
observations (2.4" xl. 9")- 

2.2.3. PdBI CO Observations 
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Figure 3. Spectral energy distribution of LABdOS focusing on 
mid-IR to mm wavelengths. New data points from this work are 
fitted with a modified blackbody SED with a dust temperature Td 
and an emissivity index /3. The error bars and upper limits are 
given at 3(7. The left inset shows the likelihood distribution of the 
fits in and /3 space. The contours represent the 1, 2 and 3 cr 
confidence intervals. We find a dust temperature, of 44, 35, 
29 K for different /3 = 1, 1.5, and 2, respectively. While there is a 
degeneracy between T^j and /3, we find that the model with colder 
dust temperature (Td ~ 30 K) provides a better fit to the data. 

Observations for SSA22-LAB01 were carried out us- 
ing the PdBI in May - June 2002 (112 GHz, tuned on 
the CO J=4-3 line at z = 3.1; Project ID: M021) and 
in June - July 2010 (84 GHz, tuned on the CO J=3-2 
line. Project ID: U045). Data were collected in com- 
pact D configuration, with baselines ranging between 24 
and 113 m. Phase calibrators were observed every 20 
minutes. The rms phase errors are 20°. Primary am- 
plitude calibrators were 3C454.3 (variable, at ~ 25 Jy at 
3mm during our observations), MWC 349 (not variable, 
^ 1 Jy) and 2223-052 (~ 3.5 Jy, also used as phase cali- 
brator). Data at 112 GHz (84 GHz) were collected using 
the narrow band (WideX) receiver. System temperatures 
ranged between 120 and 250 K. Typical uncertainties in 
the flux scales and overall calibration are about 10%. The 
data processing program used water vapor monitoring re- 
ceivers at 22 GHz on each antenna to correct amplitudes 
and phases for short-term changes in atmospheric water 
vapor. We weighted visibilities according to the inverse 
square of the system temperature, and natural weights 
were applied when creating the maps, in order to boost 
the sensitivity. The primary beam size of PdBI is 56" 
(42") at 84 GHz (112 GHz). No cleaning was applied. In 
the case of WideX observation with a large bandwidth, 
we also place an upper limit on the continuum (crms = 
0.15 mJy (3cr). We summarize the CO observations in 
Table H 

3. RESULTS 
3.1. Spectral Energy Distribution of LABdOS 

Using our new FIR photometry from Spitzer MIPS, 
Herschel SPIRE and MAMBO-2, we compare the spec- 
tral energy distribution (SED) of LABd05 with those of 
starbursts and AGNs in order to measure the bolometric 
luminosity and investigate whether we can discriminate 
possible energy sources of the blob. 

Figure [2] shows the SED of LABd05 from optical to 
radio wavelength in the observed frame. The squares 
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and large c ircles represent the flux measurements from 
iDev et all (|2005f ) and new constraints from this study, 
respectively. We also show the Spitzer Infrared Spec- 
trograph (IRS ) spec trum (Aobs = 19.5/im- 38/im) from 
iColbert et al.l (j2011| ). which reveals that strong poly- 
cyclic aromatic hydrocarbon (PAH) features contribute 
^50% of MIPS 24/im flux. For comparison, we show 
the SED templates of IRAS12112 +0305 and the Seyfert 
1 galaxy Mrk 231 that are redshifted to z = 2.66 and 
scaled to match the Infrared Array Camera (IRAC) 
3.6/im, 4.5/im, and 5.8/im flux densities (rest-frame 
l^m-l.S^fm). IRAS12112+0305 and Mrk 231 rep- 
resent local analogs whose IR SEDs are dominated 
by star-formation and AGNs, respectively. We select 
IRAS121 12 4-0305 as a representativ e template from 11 
LIRG/ULIRG templates compiled bv lRieke et aD (|2009l ) 
by searching for a template that best matches our data. 
The Mrk 231 template is also selected by searching for 
the best-fit templ ates from the AGN- dominated tem- 
plates compiled bv lPolletta et al.l (|2007[ ). 

We find that the AGN template (Mrk 231) is able to 
fit the full SED of LABd05 from the rest-frame l/im to 
350^m reasonably well, but additional PAH emission is 
required to reproduce the strong 7.7^m features in the 
observed Spitzer IRS spectrum, and the Mrk 231 tem- 
plate slightly under-predicts the long-wavelength FIR 
fiux densities. On the other hand, the starburst tem- 
plate overestimates the FIR luminosity by a factor of « 
5 and fails to fit the I RAC S.Omui and M IPS 24/im data 
points as discussed by IDev et al.l (j2005D . Clearly, more 
sophisticated SED modeling is required and observations 
around the rest-frame 2/im-10//m are essential for fully 
constraining the SED. 

Since our new (sub)mm observations cover both sides 
of the thermal dust peak, wc can constrain the dust tem- 
perature and the FIR luminosity of L ABdOS. We fit the 
data with a modified blackbody SED (lHildebrandlll983D 
with the following functional form at hv/kTd ^ 1: 



Observed wavelength [/^m] 
10000 1000 



100 



oc 



exp {hv/kTd) — 1 ' 



(1) 



where Td is the effective dust temperature, (3 is the dust 
emissivity index with 1 ^ ^ 2. In Figure [31 we show 
the best-fit models and the likelihood distribution of Td 
and /3 parameters. While Td and f3 are degenerate, we 
find that the model with the higher /3 and lower dust 
temperature Td shows better agreement with the data. 
Insufficient data on the shorter wavelength side of the 
peak prevents us from constraining the lower/upper lim- 
its on Td/P, so we choose /? = 2 as a nominal value and 
determine the best-fit dust temperature Td = 29^^ K. 
For reference, we also show the fits for fixed /3 = 1 and 
1.5 which results in Td ~ 44^3 and 35^2 K, respectively. 

Using these SED fits, we derive the far-IR and bolo- 
metric luminosity (Lfir and iboi) of LABdOS. Wc find 
iFiR(40 - lOOO/tm) ^ (4.0 ±0.52) xIO^^^q and Lboi = 
(8.6 ± 1.1) X 10^2 Lq consistent with SMGs and ultra- 
luminous infrared galaxies (ULIRGs). To obtain the 
bolometric luminosity, we apply the Lboi ~ 2.15 Lfir 
correction based on the Mrk 231 SED template. Note 
that the SPIRE 350Aim flux density, vL^ = (3.5 ±0.6) 
xlO^^L0, alone can roughly constrain the Lfir since 
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Figure 4. FIR flux densities of SSA22-LAB01. Different symbols 
represent the (sub) mm observations with different instruments and 
beams sizes (see Table O- The filled and open symbols indicate 
the flux measurements from single-dish and interferometer obser- 
vations, respectively. All the uncertainties and upper limits are in 
3(7 level. For illustration purposes, the solid and dot-dashed lines 
represent the modified blackbody SEDs normalized at the SssOfjm 
= 6mJy with (T^j, 13) = (40K, 2) and (30 K, 1.5), respectively. The 
upper limits from all other instruments disagree with the SCUBA 
850/im measurement. 

the SPIRE band corresponds to the peak of the spectral 
energy distribution. We estimate the dust mass, Md, 
using 

_ 1 S^Dl 

where Si, is the observed flux density, is the lumi- 
nosity distance, n'jf^^ is the rest frame dust mass absorp- 
tion coefficient, and B{i', Td) is the Plank function at the 
rest frame. We adopt Td — 29 K and (3 — 2 from our 
best-fit model and Ki25fj,m = 26.4 cm^ g~^ ([Dunne et al.l 
I2003D . which corresponds to Kgso^m = 0.57 cm^ g~^ as- 
suming K (X v^ . We find a total dust mass of Mjust = 
(1.4 ± 0.7) X 10^ Mq. Note that this dust mass estimate 
strongly depends on the choice of k™'^' and T^. The ab- 
sorption coefficient Kgso^ni is uncertain to a factor of ~2; 
e.g., Kssopm = 0.77 cm^ g-i ([Dunne et al1[2000l [20Tll) or 
Ks.-sn/t TTi = 0.38 cm^ g~^ (Milky Way dust model: [Praind 
I2003f ). Or if we adopt the lower temperature = 35K 
for (3 = 1.5, the dust mass will decrease by a factor of 
2.2. 

3.2. Non-detection of Continuum Emission in 
SSA22-LAB01 

SSA22-LAB01 was not detected down to /SsTO/im = 
8.1 - 12mJy beam"! (3cr) with new APEX/LABOCA 
single-dish observations which have similar beam size 
(19") t o the previous SCUB A 850//m observations 
(15"). [Chapman eFall ([2004D reported that SSA22- 
LABOl contains a bright submm galaxy: SMM 
J221725.97+001238.9 with the flux densities of S'sso^tm 
= 16.8 ± 2.9 mJy and S'45o^m = 45.1 ± 15.1 mJy, respec- 
tivelyS If we assume simple model SEDs with (T^, /3) 
= (40K, 2) and (30K, 1.5), we expect an 870/xm fiux 
density of S'sro/^m = 15.7 - 16.1 mJy, which would be 
detected by LABOCA with 4 -6a significance depend- 
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Figure 5. {Left) PdBI 1.2 mm continuum map of SSA22-LAB01. The syntliesized beam has a FWHM of 2'.'4xl'.'6 {P.A.=14.9°). The 
phase center of the PdBI is marked with a large cross — its size indicates the positional uncertainty of the LABOCA pointing. The 
LABOCA beam (F WHM = 19") is represented by the dashed circle. The small crosses represent the locations of the Spitzer IRAC sources 
IjGeach et al.|[200W ) . No significant continuum source is detected above Scr level at the location of these galaxies. (Middle an d right) PdBI 
1.2 mm continuum contours superimposed on the Subaru/Suprimo-Cam hya and BV broadband images I IMatsuda et al.ll2004l) . respectively. 
The contours show the —3, —2, —1, 1, 2, Scr with a = 0.15 mjy beam~^. 



ing on the LABOCA error estimates fi )2.2.ip . This non- 
detection is consistent with essentially all measurements 
of the source in the literature, in pa rticular the SMA 
observations by IMatsuda et al.l (|2007D . In that study, 
SSA22-LAB01 was also not detected down to S'sso/xm 
= 4.2 mJy beam^^ (3cr) with a ~2" spatial resolution. 
They proposed, as explanation, that the dust emission 
(reported by SCUBA) is extended on spatial scales larger 
than 4" such that the interferometer observation would 
resolve out this smooth component. However, the origi- 
nal SCUBA observation is not verified by our LABOCA 
observation, for which this argument does not hold. Fur- 
thermore, an ASTE-AzTEC single-dish observation with 
a larger beam size (FWHM w 28") than SCUBA (15") 
also failed to d etect dust continuum down to mJ y 
(3f7) at 1.1mm (|Kohno et al.l[200l[Tamira et al.l[2009l ). 
thus contradicting the SCUBA measurement at the 4.5cr 
confidence level. In Figure |31 we show the upper lim- 
its derived from our LABOCA and PdBI observations 
with the existing (sub)mm observations of the SSA22- 
LABOl. These observations are summarized in Tabled 
The solid and dot-dashed lines represent the two model 
SEDs normalized at S'sso^m = 6mJy, showing the pos- 
sible range of the extrapolated flux densities at various 
wavelength. The upper limits from all other (sub)mm 
observations (both single-dish and interferometer) are in- 
consistent with the SCUBA measurements. Therefore, 
we conclude that the previous SCUBA measurements are 
not reliable. 

Given that SSA22-LAB01 contains multiple galaxies 
and that the total star formation rate within the blob is 
at least '-22OM0yr-^ (iJylatsuda et al. 2003), it is criti- 
cal to identify galaxies by their dust emission within the 
Lya halo which produce most of bolometric luminosity. 
In order to detect dust emission from individual galax- 

^ Slightly different values for the SCUBA 850itm a nd 450/im 
fluxes are reported in lChapman et al.l l l200ll. [200l. 120051 ') . Since all 
the flux measurements agree wi th each other within u ncertainties, 
wo adopt the values reported bv lChapman et al.l l|2004l ) throughout 
the paper. 



ies within SSA22-LAB01, we compare our deep PdBI 
1.2 mm continuum map with the locat ions of the UV 
and S pitze r IRAC sources ide ntified by iChapman et al.l 
(|2004l ) and lGeach et al.l ()2007f ). In Figure El we show the 
PdBI 1.2 mm continuum map of SSA22-LAB01 and over- 
lay the conto urs of this map on th e Lya and BV broad- 
band images ([Matsuda et al.ll2004[ ). The pointing center 
of the LABOCA observation is marked with a large cross 
(positional uncertainty of 3'.'6) and a dashed circle gives 
the beam size of 19". This pointing center is the same 
as the phase center of the PdBI observations. We show 
location s of the Spitzer/IR AC sources marked as LAB 
1 abcde (jCeach et all 120071 ) in Figure [H Out of these 
five IRAC sources, three sources (LAB 1 a, b, and e) are 
likely associated with the SSA22-LAB01 (jOeach et al.l 
120071: lUchimoto et al.ir2008f ). We find that no significant 
emission is detected from these three galaxies while there 
is a -4ct peak (R.A. = 22in7"26!06, deck = 00°12'35"1) 
between LAB la and 16. Furthermore, none of ~3(t 
peaks within the FOV of our 1.2 mm map agrees with 
other UV sources in the field (right panel in Figure [5]) . 
Therefore, we conclude that none of the known optical or 
near-IR emitting galaxies within the SSA22-LAB01 Lya 
halo are detected at 1.2 mm down to 0.45 mJy beam~^ 
(3cr). Future observations with the Atacama Large Mil- 
limeter/submillimeter Array (ALMA) will be required to 
verify this tentative ■~4(t peak. 

Without any formal detection at (sub) mm wave- 
lengths, the FIR luminosity of SSA22-LAB01 remains 
highly uncertain. We place upper limits on the FIR lu- 
minosity of SSA22-LAB01 using simple model SEDs to 
provide the possible range of Lfir- If we normalize the 
SEDs at the 870/im 3a upper limit, we obtain ipiR < 
3.1xl0^^io and 1.2x10^^^© depending on the choice 
of the model SED, (T^, /?) = (30 K, 1.5) and (40 K, 2), 
respectively. If we adopt our 1.2 mm PdBI measurement, 
the individual galaxies within the blob have Lfir < (0.39 
- 2.1) xlO^^ (3cr). 

In addition to the PdBI 1.2 mm continuum observa- 
tion, we also place an upper limit on the 3 mm flux den- 
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Figure 6. {Top) IRAM-30m CO spectra for CO J=3-2 and CO J=5-4 transitions for LABdOS. Both spectra have velocity resolutions 
of 5v = 40kms~^. The thick and thin horizontal bars represent the expected velocity range estimated from the Lyo and He II emission 
line in the optical l|Dev et al.ll2005l ). {Bottom) PdBI CO spectra for CO J =3-2 and CO J = 4-3 transitions of SSA22-LAB01 with 5v = 
lOOkms"^ and 27kms~^, respectively. Note that the scale s of y-axis are diffe rent between two panels. The thick horizontal bars represent 
the expected velocity r ange estimated from t he hya line l|Bower et al.|[200i ). The vertical arrows indicate the location of the previous 
tentative CO detection ijChapman et al1l2004l 1. No CO line is detected in either case. 

sity using the PdBI CO observation (see Section 1331 and 
Table [3]). No continuum source is detected in the inte- 
grated channel map over 4 GHz bandwidth, resulting in 
3(7 limit of 0.15 mJy beam^^ with a beam size of 8'.'4 x 
5'.'4. 



3.3. Molecular Gas Content of LABdOS and 
SSA22-LAB01 

In Figure [H we show the spectra of the CO lines for 
the LABd05 and the SSA22-LAB01 obtained from the 
IRAM-30m and PdBI, respectively For LABdOS, we 
show the spectra for CO J =3-2 and CO J =5-4 lines 
with a velocity resolution of 5v = 40kms~^. For SSA22- 
LABOl, we show CO J=3-2 and CO J=4-3 transitions 
with 5v = lOOkms"^ and 27kms~^, respectively. The 
rms noise per channel ranges from 0.5 to 2.8 mJy beam~^ 
(Tabic [21) depending on the instruments and the transi- 
tions. 

Because of the width of spectral bandwidth 
(^2000 km s^^), knowing the accurate redshift is 
critical in this CO detection experiment. In Figure |6l 
we show the velocity range of CO lines that corresponds 
to the range spanne d by Lya: z = 2.656 ± 0.006 for 
LABd05 (Dcv et all [20051) and z = 3.102 ± 0.005 fo r 
SSA22-LAB01 (jBower et al.l[200l IMatsuda et al.l[200l . 
Given that Lya emission lines from star-forming galaxies 
at z = 2 — 3 could be redshifted against the optically 
thin nebular lines (i.e., systemic velocity) by 250 - 
lOOOkms-i fe.g.. lPettini"eraI|[200lt ISteidel et al.][200l 



l2010f l due to either galactic-scale outflows or absorption 
by the intervening intergalactic medium, we note that 
CO lines, which are a,lso expected to emit near the 
systemic velocity (e.g.. iGreve et al.ll2005D . could be lo- 
cated blueward (higher I'ohs) of the horizontal bars. For 
LABd05, we also show the velocity range corresponding 
to the optically thin He II A 1640 emission that originates 
from the central region of the Lya nebula. At least 
for LABd05, these two estimates for CO line centers 
are in good agreement, therefore wc conclude that it is 
unlikely that any potential CO line would fall outside 
of the spectral bandwidth covered by our IRAM-30m 
observations. We note that, for SSA22-LAB01, the 
CO J=4-3 spectra obtained with the old narrow-band 
receivers at PdBI were taken with a limited bandwidth 
and no other optically thin emission lines were available 
to double-check the systemic redshift. On the other 
hand, the PdBI CO J=3-2 spectrum obtained with the 
new WideX correlator has a very wide spectral coverage 
(4GHz; 14000kms~^), enough to cover any plausible 
velocity offsets between Lya and CO lines. We conclude 
that no significant CO emission is detected above the 
rms noise near the expected velocities in both systems. 

In the case of SSA22-LAB01, we also inspect whether 
any tentative source can be detected in the integrated 
channel maps. In Figure jT] wc show the channel maps 
summed over AV = 1000 km s~^ intervals centered at 
-1500, -500, 500, 1500kms"i from the Lya velocity 
center. While there could be a possible 3cr detection near 
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the phase center in the v — — SOOkms"^ channel map 
(upper right panel of Figure ^ , this marginal detection 
in the integrated map corresponds to only ^1.7 a signal 
per lOOkms"^ channels. Because we are not able to 
extract a reliable spectrum from this tentative source, we 
conclude that no significant CO line is detected above 
the 3 X rms noise near the expected velocities. More 
sensitive observations using ALMA will test whether or 
not this tentative emission is real. 

Our new PdBI observation for SSA22-LAB01 rules out 
the previous ten t ative CO detection at high significance. 
IChapman et al.l (|2004[ ) reported a tentative (3.2 cr) de- 
tection of a CO J=4-3 line from SSA22-LAB01 using 
the Owens Valley Radio Observatory Millimeter Array. 
The reported intensity is Si, ^ 10 mJy at the peak with 
a line width of ^400kms~^. In Figure IHl we indicate the 
location of earlier tentative CO detection with arrows. 
As described in detail below, our 3ct limit for the inte- 
grated CO J=4-3 flux from PdBI is AV < 0.62 Jy 
kms~^, i.e., we exclude the previous detection at high 
significance (12(t). 

Using these non-detections, we put constraints on CO 
line lumi nosities and molecular gas mass of the blobs. 
Following ISolomon fc Vanden BoutI (|2005[ ). the CO lu- 
minosity L'qq (in Kkms~^pc^) is given by 

L'co = 3.25 X 10^ Sco AV v'^^ Dl (1 + z)-^ , (3) 

where fobs is the observing frequency (in GHz), is 
the luminosity distance (in Mpc) to the source at a red- 
shift z. For our non-detection, we adopt a 3ct upper limit 
on the velocity-integrated flux S'coAT^ = 3(AT//(St>)i/2 
[5v (Trnis) in Jykms"^, where AV is the CO line width, 
where 5v is the channel bandwidth, and u-cms is the rms 
noise value per channel, respectively. As the CO line 
width is unknown, we adopt AV = 400kms~^, consis- 
tent with FWHM = 300 - 450 km s^^ measured from the 
optically thin Ha lines in two other Lyg blobs in th e 
Extended Chandra Deep Field South ()Yang et al.ll2011[ ). 
For LABd05, we obtain 3cr CO line luminosity limits, 
ico(5-4) < 1.37X 10^°(AT//400)i/2 Kkms'^pc^ and 
ico(3-2) < 3.90 X 101° (Ay/400)i/2 Kkms'^pc^. In 
the case of the SSA22-LAB01, we obtain the upper lim- 
its of Lco(4-3) < 1-6 X 10i°(Ay/400)i/2 and 2.co(3-2) 
< 1.5 X 10i°(AF/400)i/2 Kkms-ipc^. 

The upper limits for the total molecular gas mass, 
M(H2), can be derived using a CO-to-H2 conversion 
factor between CO J=1~0 luminosity and molecular 
gas mass, Xcq. Here we adopt Xqq ~ 0.8 Mq 
(Kkms~i pc^)~i, which was proposed as appropriate for 
starburst environments found in ul tra-luminous infrared 
galaxies (jPownes fc Solomonl[T998l) . If we assume a con- 
stant brightness temperature for the different CO tran- 
sitions, i.e., CO J= 1-0 luminosity as defined by Eq. ([3)) 
is the same as those of upper (J+l) — > J transitions, the 
upper limits on the molecular gas mass can be given as 
M(H2) < 3.1 and 1.2 x 10^'^Mq{AV/A00)^^^{Xco/0.8) 
for LABd05 and SSA22-LAB01, respectively. Note that, 
for submm galaxies and QSOs in a similar rcdshift range, 
the gas has been found to be thermalized at least up to 
J=3-2 transition as the r atio i^co(3-2)/-^co (i-o) is close 
to unity (w 0.8) (e.g., IWeiss et al.ll2007t) . Therefore, 
we use only the J=3-2 transition to estimate limits on 




Figure 7. Integrated channel map of the CO J =3— 2 line for 
SSA22-LAB01. Each channel map is integrated over a 1000 km 
bandwidth centered at the velocities shown in the upper left cor- 
ners. The contours are at —3, —2, —1, 1, 2, and Scrims, where 
(Trms is the rms noise per lOOOkms-l velocity width (0.172 mJy 
beam~^). The phase center is marked with a cross in each panel. 
While there is a possible Scr detection near the phase center at 
the V = —500 km s^^ channel, no significant CO line is detected at 
above So" level per lOOkms"^ channel in that bin. 

Lco(i-o) and M(H2). Wc note that our final A/(H2) es- 
timates are strongly dependent not only on the intrinsic 
excitation, but also on the choice of Xqo- For exam- 
ple, if wc adopted Xco ~ 4.5 Mq (Kkms^^ pc^)"^, the 
Milky Way value, the resulting molecular gas mass limits 
would increase by a factor of ~^5. 

3.4. LpiR-L'co Correlation 

With the FIR luminosity and the upper limit on CO 
luminosity in hand, we investigate whether Leo a-nd 
Lpm of the LABd05 are consistent with the known 
ipiR-ico scaling relations. For SSA22-LAB01, the 
ratio between LpiR and Leo remains unconstrained 
so far. For nearby starburst and spiral galaxies with 
LpiR ^ IO^^Loj the correlation has a form, logLpiR 
= (1 .26 ± 0.08) X logL^o - 0-81 dCao fc SolomonI 
[200l . When ULIRGs, SMGs, radio galaxies (RGs), and 
QSOs are included, the relation becomes slightly steeper, 
logLpiR = (1.39 ± 0.05) X logL^o - 1-76, but still 
holds over 5 orders of magnitude in LpiR (jRiechers et al.l 
[20M ). For LABd05 with LpiR = (4.0 ±0.5) x lO^^ Lq, 
we expect L'co = (2.2±1.9) x 10^" Kkms"! pc^ for the 
ULIRG/SMG/RG/QSO relation. Therefore, we find 
that our 3(7 upper limit on L'qq for LABd05 is consistent 
with the known ipiR-L'co scaling relations. In terms 
of the continuum-to-line luminosity ratio, we find that 
Lpm/L'co > 100, which is consistent with the LpiR/ico 
= 125-240 L0(Kkms~^ pc^)~^ for the galaxies with 
LpiR ^ 10^^^^^ Lq. Note that due to the non- linear re- 
lation between LpiR and L'cOi the ratio ipiR/ico has 
a large spread depending on ipiR. 
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4. SUMMARY AND CONCLUDING REMARKS 

We have obtained IR and (sub)mni observations of 
the two best-studied Lya blobs, LABd05 and SSA22- 
LABOl; in order to constrain their energy budgets (i.e., 
the bolometric luminosity) and dust properties. We 
find that LABd05 has a high FIR luminosity, Lfir= 
4.0 X 10^^ Lq, comparable to values found for high-z 
SMGs and ULIRGs. The NIR-to-FIR SED of LABdOS 
is well described by the AGN-starburst composite tem- 
plate (Mrk 231). For SSA22-LAB01, no STO^im contin- 
uum is detected down to 8.1-12 mJy beam~^ (Scr) in the 
LABOCA single-dish observation contrary to the origi- 
nally reported SCUBA measurement. To detect the dust 
emission from individual galaxies within SSA22-LAB01, 
we obtained a very deep 1.2 mm observation of SSA22- 
LABOl. No 1.2mm continuum is detected down to w 
0.45 mJy beam~^ {3a) at ~2" resolution. Combined 
with the existing extensive (sub)mm observations from 
literature, we con clude that the previo usly published 
SCUBA detection (iChapman et aLlboOl) is not rehable. 
We place 3cr upper limits of Lfir < (3.1 - 12) x 10^^ 
Lq for the total FIR luminosity of the blob, and Lfir < 
(0.39 - 2.1) X 10^^ Lq for the individual galaxies within 
the blob. 

To investigate the molecular gas content in Lya blobs, 
we carried out a sensitive search for CO lines in the two 
Lya blobs. No CO line is detected down to an inte- 
grated flux limit of S^AV < 0.25 - 1.0 Jykms"! con- 
straining the molecular mass to be less than M(H2) « 
(1-3) X IO^^Mq (3cr limit), assuming a constant bright- 
ness temperature and a CO-to-H2 conversion factor for 
the starburst galaxies. The non-detections from our 
new search exclude the previous tentative (~3cr) de- 
tection of th e CO J =4-3 l i ne tow ard SSA22-LAB01 
reported by IChapman et al.l (|2004[ ) with a high sig- 
nificance (12cr). We find that the FIR-to-CO lumi- 
nosity ratio of Lpm/L'co ^ 100 L0(K kms~^ pc^)~^ 
for LABd05 is consistent with the scaling relations 
of ULIRG/SMG/QSO/RG. While our sensitive CO 
searches already place interesting limits on the bright- 
est Lya blobs, future observations with ALMA will al- 
low us to routinely detect CO in similar systems (0.1 
Jykms"^ for ■^l hour integration) beyond the sensitiv- 
ity limits that are accessible today. 
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Table 1 

Photometry of LABdOS 



Facility 


-^obs 




Note 




{fim) 


(mjy) 




Spitzer/MIPS 


24 


0.856 ± 0.005 


Dcy et al. (2005) 




70 


< 9 


this study 




160 


< 51 


this study 


CSO SHARC-II 


350 


37 ± 13 


Bussmann et al. (2009) 


Herschel SPIRE 


250 


18.8 ± 5.2 


this study 




350 


26.9 ± 5.1 


this study 


IRAM MAMBO-2 


1200 


2.76 ± 0.35 


this study 



^ 3(T upper limits if non-dctcction. 



Table 2 

CO Line Observations for LABd05 and SSA22-LAB01 



Source 




Transition 


I'obs 




5v^ 


L'co {j+i^jr 


M(H2)<^ 








(GHz) 


(mJy) 


(kms-i) 


(10^° Kkms-ipc2) 


(10i« Mq) 


LABd05 


2.656 ±0.006 


3-2 


94.614 


2.8 


42.0 


< 3.90 


< 3.12 






5-4 


157.675 


2.8 


42.0 


< 1.37 




SSA22-LAB01 


3.102 ±0.005 


3-2 


84.299 


0.53 


100 


< 1.47 


< 1.17 






4-3 


112.300 


2.0 


26.7 


< 1.64 





" Redshifts determined from Lya lines <Dev et al JI2005I : IBower et al.1120041 '1 . 
^ Channel bandwidth. 

^ 3cr upper limits assuming CO linewidth AV — 400 km s^^. 

3(7 upper limits assuming Xqo — 0.8 Mq (Kkms""^ pc^)~^ and constant brightness temperatures. 



Table 3 

FIR Observations of SSA22-LAB01 



Facility Wavelength Beam size S^^ Reference'' 

(^m) (arcsec) (™Jy) 

SCUBA 850 15 16.8 ± 2.9 (1) 
450 8 45.1 ± 15.1 

AzTEC 1100 28 < 3.0 (2) 
SMA 880 2.5x1.9 < 4.2 (3) 

LABOCA 870 19 < 12.0 this study 

PdBI 1250 2.4x1.6 < 0.45 this study 

PdBI 3500 8.4x5.4 < 0.15 this study 

^ 3(7 upper limits if non-detection. 

^ (1) [Chapman et al] I I2001I, [200l. 120051), (2 ) IKohno et al] 1 120081') : 
ITamura et al.l 120091) , (3) IMatsuda et al.l I I2007I ) 
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